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The role of compressibility in solar wind plasma turbulence. 
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Incompressible Magneto-hydrodynamics is often assumed to describe solar wind turbulence. We 
use extended self similarity to reveal scaling in structure functions of density fluctuations in the 
solar wind. Obtained scaling is then compared with that found in the inertial range of quantities 
identified as passive scalars in other turbulent systems. We find that these are not coincident. This 
implies that either solar wind turbulence is compressible, or that straightforward comparison of 
structure functions does not adequately capture its inertial range properties. 

PACS numbers: Valid PACS appear here 



The supersonic and super-Alfvenic flow of the solar 
wind offers a unique opportunity to investigate the prop- 
erties of the magnetized and turbulent plasma. The tran- 
sition from a laminar to turbulent flow requires large 
Reynolds number Re = LU jv, and its magnetic coun- 
terpart R m = LU/r], where L is the energy injection 
scale length, U is the velocity difference on scale L, v 
is the viscosity and 77 is the magnetic diffusivity. Es- 
timations of the hydrodynamic and magnetic Reynolds 
numbers in the solar wind exceed 10 8 [l], Q compared 
with Re w Rm « 10 4 achieved in direct numerical sim- 
ulations (DNS) [J and just few hundreds in some 
magnetized liquid laboratory experiments 0, Q- The 
presence of turbulence in the solar wind is strongly sug- 
gested by numerous observations. These include power 
law power spectra with —5/3 Kolmogorov-like slopes in 
the kinetic and magnetic energy densities (eg., H,S0) 
and non-Gaussian Probability Density Functions (PDFs) 
(e.g., jTTL IT2I IT3I Ilil] ) found for fluctuations in the veloc- 
ity and the magnetic field. 

These observations imply that solar wind turbulence 
shares many of its statistical properties with incompress- 
ible isotropic hydrodynamic turbulence 0,0>EJ As a 
result, the turbulent dynamics of the solar wind is often 
modelled assuming that the plasma density is constant. 
This assumption of incompressibility is particularly con- 
venient in analytical and numerical studies of magneto- 
hydrodynamic (MHD) turbulence. The assumption of 
incompressibility also appears to be in good agreement 
with the results of compressive MHD simulations where 
generation of compressive modes from Alfvenic turbu- 
lence was found to be suppressed 2} . In the context of 
the solar wind, incompressibility has been suggested to 
be a reasonable approximation for plasma in fast wind 
streams pH, 111. Il9) . Considerable progress has been made 
by treating the solar wind as dominated by Alfvenic fluc- 
tuations (e.g., 113, El H3, and with nearly incom- 
pressible magnetohydrodynamic theory j2(|. 

However, statistical features associated with turbu- 
lence have also been identified in the density fluctuations 
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FIG. 1: Conditioned structure functions of fluctuations in 
the density p for slow solar wind. 



derived from the solar wind observations. A k~ 5 / 3 scal- 
ing in the omnidirectional wave number spectrum of elec- 
tron and proton density has been reported as early as the 
1970s (see [24| and the references therein). Models of the 
solar wind can offer possible mechanisms for generating 
strong density fluctuations close to the heliospheric cur- 
rent sheets [23. Recent studies suggest that, at least in 
the slow solar wind, the coupling of Alfvenic and mag- 
netosonic MHD modes could not fully describe the na- 
ture of these fluctuations. For example, the ratio Sp/ (p) 
was found to be nearly constant and independent from 
the amplitude of the magnetic field fluctuations SB /(B) 
in the slow solar wind[25|]. More recently the density 
fluctuations were also found to exhibit non trivial, ap- 
proximately self-affine scaling^4| similar to that found 
in other solar wind bulk plasma parameters. 

In this Letter we examine the scaling properties of the 
proton density p in fast and in slow solar wind. We focus 
on a comparison of the scaling properties of fluctuations 
in density with that of various passive scalars identified 
in experiments, and direct numerical simulations, of fluid 
turbulence. Previously, these passive scalars have been 
shown to exhibit scaling that corresponds closely to that 
of the magnitude of the magnetic field B in the solar 
wind, although the data interval was not ordered by solar 
wind speed 28] . By assuming incompressibility it can be 
argued [2|| that B should act as a passive scalar. We 
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FIG. 2: ESS of conditioned structure functions of fluctuations 
in the density p for slow solar wind. 
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FIG. 3: ESS of conditioned structure functions of fluctuations 
in the density p for fast solar wind. 

will for completeness also examine the scaling properties 
of magnetic field magnitude in fast and slow solar wind. 

It is instructive to recall that the dynamics of a passive 
scalar T = T(x, t) in a velocity field v(x, t) is given by 
the advection equation 

d t T = -(v V)T + kV 2 T. (1) 

where n is the diffusivity. Dynamical and statistical prop- 
erties of a passive scalar are more tractable analytically, 
as compared to active fields like velocity, since the equa- 
tion Q is linear in T . In the case of hydrodynamics, 
theoretical predictions have also been verified, to some 
extent, experimentally using tracer particles that do not 
disturb the flow [2^. Let us consider the compressible 
MHD equations for the evolution of the magnetic field 
and the density: 

d t B = Vx(vxB)+??V 2 B and %»+V-(v/?) = 0. (2) 

Given the assumption of incompressibility (V • v = 0), 
it immediately follows that p should behave as a passive 
scalar in the turbulent solar wind flow: 

d t p=-(vW)p. (3) 

To investigate the scaling properties of the density p 
and the magnetic field magnitude B we use 64 seconds 
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FIG. 4: ESS of conditioned structure functions of fluctuations 
in the magnetic field magnitude B for slow solar wind. 



averaged data from the ACE spacecraft |3Jj] set spanning 
from 01/01/1998 to 12/31/2001. This interval includes 
dates previously considered in Ref. 28] . The slow and fast 
solar wind are known to exhibit distinct phenomenology 
(e.g., We thus split the data into slow and fast 

solar wind sets using 450km/s wind speed as a separation 
criteria. The resulting data sets consist of ~ 1 x 10 6 sam- 
ples for the slow wind and ~ 0.6 x 10 6 for the fast wind. 
We apply structure function|3l| analysis to fluctuations 
in density for slow and fast wind separately. Generalized 
structure functions S m of fluctuations in say, the density 
p(t), on timescale r are defined as moments m through 
S m (r) = (\p(t + r) — p{t)\™) where the ensemble (. . .} is 
taken in the time domain [33] ■ If scaling is present in the 
time series we expect these to show a power law depen- 
dence on the temporal scale r, i.e., S m oc r^ m \ 

Finite, experimental data sets include a small number 
of extreme events (outliers) that, due to poor statistics, 
may obscure the correct scaling of the high order mo- 
ments. Here, we will exclude these events by the use 
of conditioning [17], This approach puts a limit on the 
range of fluctuations used in computing structure func- 
tions. This limit is varied with the temporal scale r to 
account for the growth of range with temporal scale in the 
signal. In our case we defined this threshold as 15ct(t), 
where <t(t) is a standard deviation of fluctuations on tem- 
poral scale t. We stress that conditioning improves the 
scaling where it already exists but does not enforce it 
on the investigated data, if the applied threshold is suf- 
ficiently large. In practice, for the limit chosen here, we 
eliminate less then 1% of the data points. 

Figure^Jshows the structure functions S m plotted ver- 
sus t on logarithmic axes for orders 1 < m < 4 for fluctu- 
ations in density in the slow solar wind. The plot shows 
a scaling region extending from r ~ 10 minutes to r ~ 3 
hours (~ 1.5 decades on the logarithmic scale). The qual- 
ity of the scaling deteriorates when fast wind streams are 
considered. The scaling regions can be still identified 
but they extend only < 1 decade from ~ 10 minutes to 
1 hour. 

It has been empirically shown that Extended Self Sim- 
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FIG. 5: ESS of conditioned structure functions of fluctuations 
in the magnetic field magnitude B for fast solar wind. 
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FIG. 6: Scaling exponents of the conditioned structure 
functions for SB (squares) and Sp (circles) in slow solar 
wind. Dashed line and diamonds correspond to DNS|3| and 
experimental [35il results, respectively. 



ilarity (ESS) approach can consider ably extend the re- 
gion of scaling in structure functions 33] . The method is 

based on the assumption that S rn (T) oc Sp^(r). This 
suggests that scaling should emerge when the quantity 
S m is plotted as a function of S p . The scaling exponent 
£ can then be obtained straightforwardly from the rela- 
tion C(m) = Q{p)fl{m)- For the density fluctuations in 
the solar wind, we find that S3 is close to unity and we 
plot S m versus S3 on logarithmic axes for fluctuations in 
the density in slow and fast solar wind, in figures |2] and |3] 
respectively. These figures demonstrate that we can ex- 
tend scaling in the density in both the slow and the fast 
wind streams to over 2 decades when ESS is applied. 

ESS was previously applied to the magnitude of mag- 
netic field in an undifferentiated interval of solar wind|28| 
and for completeness we give the analysis in slow and fast 
wind here. Figures 0] and |S] show S m versus 53 on log- 
arithmic axes for fluctuations in the magnitude of mag- 
netic field. Previously, the scaling exponents were ob- 
tained by considering S m versus 64 |28( and we have ver- 
ified that these closely correspond to the values found 
here. The local slopes of the third order structure func- 
tions 5*3 then yield an estimate of the exponent £(3) and 
these are detailed in Table U 

We can now directly compare the scaling found for 
the density p with that identified for quantities acting 
as passive scalars and that of the magnetic field mag- 
nitude B. The resulting functional form of the scaling 
exponents C( TO ) f° r the density and magnetic field mag- 
nitude in slow and fast solar wind are shown in figures 
[5] and [7| respectively. For comparison, the scaling expo- 
nents obtained for passive scalars from the DNS |3J| and 
the the wind tunnel experiment |35| are also shown. We 
immediately see that whereas the exponents for magnetic 
field magnitude and the passive scalars fall close to each 
other on these plots (as also reported in j2^) they are 
distinct from those obtained for the solar wind density in 
both slow and fast solar wind. The fluctuations in 6B, for 
slow and fast wind, and that of the passive scalars exhibit 
multi-fractal scaling. Intriguingly, the density fluctua- 
tions are nearly self-affine with scaling exponent a that 
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FIG. 7: Scaling exponents of the conditioned structure func- 
tions for SB (squares) and Sp (circles) in fast solar wind. 
Other symbols same as in figure [(J 

differs in slow and fast wind. The values of these expo- 
nents are: a p slow = 0.39 ± 0.03 and a p fast = 0.33 ± 0.03. 

This suggests one of two possible conclusions. The first 
of these is that the turbulent solar wind is compressible 
and equation does not hold, so that the density is an 
active quantity. This is rather surprising and calls into 
question the significant body of work on MHD turbu- 
lence in the solar wind that relies on the assumption of 
incompressibility (for example, l2 ^|2lll22L 23] ). This also 
invalidates the arguments in Ref. [28| required to cast the 
advection equation for B in the form of a passive scalar 
which require incompressibility, yielding: 

d t B = -(v- V)B + r/V 2 S + AS, (4) 

The second possibility is that commonality of scal- 
ing (as measured by structure functions) does not im- 
ply shared phenomenology, and as a corollary, an ab- 
sence of such a commonality does not imply distinct phe- 
nomenology. The fact that several quantities share the 
same structure functions is then coincidental rather than 
expressing some universality of fluid turbulence. This im- 
plies that the structure functions of a single quantity do 



4 



TABLE I: Exponents £(3) obtained from figures^ 



Quantity Slow/Fast 


C(3) 


Ranj 


;c [min] 


S\B\ Slow 


0.85 ±0.06 


10 


- 500 


Sp Slow 


1.16 ±0.08 


10 


- 120 


S\B\ Fast 


0.79 ±0.05 


8- 


- 200 


8p Fast 


1.01 ±0.07 


6- 


- 100 



not fully capture the phenomenology of a given turbulent 
system. This has implications for analysis of turbulence 
which has been achieved through comparison between 
(multi-fractal) models of turbulence and the data, via 
the structure functions (see, for example, [12). 

A possible resolution may be found in the form of equa- 
tions © and (J3J . It is known that the molecular diffusiv- 
ity can change the inertial range scaling properties of a 
passive scalar [^ . The observed differences could then be 
due to the presence of the diffusive term in the magnetic 
field equation Q and its absence in the density equation 
This may account for the different scaling found in 
these quantities. 

Intriguingly, equation @, when written for the mo- 
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ments of fluctuations in p l has no explicit dependence on 
the order of the moment [36| 

[d t + v(x, t) ■ V ± v(x', t) • V] (Sp) n = 0, (5) 

where (Sp) n = [p(x,t) — p(x!,t)] n . In this case fluctua- 
tions in density should be simply those imposed on the 
initial condition 5p(x, 0). These may be mediated via 
large scale coherent structures (shocks and coronal mass 
ejections). This may suggest a solar origin of these fluc- 
tuations in the density. It may thus be informative to 
attempt to relate the scaling found in the density fluctu- 
ations in the solar wind, with that of the solar corona. 
The power law scaling of X-ray flux from solar flares [3^ 
is intriguing in this regard. 
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